卒業論文発表会
1
月28
日, 2009, 福井大学工学部物理工学科量子エンタングルメント の指標とその突然の消滅
物理工学科 57山下悟士 44丹尾大輔
量子エンタングルメント
古典 bit (0 , 1)
例 コインの裏と表
量子 bit (qubit) 2 つの状態 ( | 0 i , | 1 i ) しか持たない系
例 電子のスピン
(
右回り| 0 i
、 左回り| 1 i )
2準位原子(
基底状態| 0 i
、 励起状態| 1 i )
重ね合わせの原理→ | φ i = α| 0 i + β| 1 i
• qubitが2つの場合での量子系
| 00 i , | 01 i , | 10 i , | 11 i
重ね合わせの原理| φ i = α 00 ( | 00 i + α 01 | 01 i + α 10 | 10 i + α 11 | 11 i
|α 00 | 2 + |α 01 | 2 + |α 10 | 2 + |α 11 | 2 = 1 )
ある
2qubit
の状態| Φ i = √ 1
2 | 0 i a | 1 i b + √ 1
2 | 1 i a | 0 i b
確率:1
2 · · · | 0 i a → | 1 i b
確率:
1 2 · · · | 1 i a → | 0 i b
⇓
測定結果が絡み合っている
(エンタングルメント )
このようにエンタングルした状態を、最大エンタングルメント 状態又はEPR pair( Einstein・Podolsky・Rosen
対)
という。研究の目的
•
エンタングルメント の有用性の例 量子テレポーテーション•
量子テレポーテーションの忠実度とエンタングルメント の指標(Concurrence)
量子テレポーテーションの忠実度F (Fidelity)
Concurrence
の導入• 2
準位原子系におけるConcurrence
の時間変化エンタングルメント の突然死
(Entanglement Sudden Death)
• 2
準位原子系における有限温度でのConcurrence
の時間変化 有限温度でのエンタングルメント の突然死量子テレポーテーション
目的:
Alice
に与えられた任意の状態のqubit
を、遠く離れた位置にいるBob
にテレポートしたい。
1.
最大エンタングルメント 状態|Φi
abの2qubit (qubit a , qubit b)
を作り、Aliceがqubit a、Bob
がqubit b
を持ち遠く離れる。2. Alice
に、Bob
へ送りたい状態| φ i
cのqubit c
を与える。3. Alice
はqubit a
とqubit c
に対して、ある測定を行い測定結果(2bit)
をBob
に伝える。4. Bob
はAlice
からの情報に応じて、qubitb
に対して特定の変換を行い|φi
cを得る.qubit A qubit B
Alice Bob
エンタングルメント
qubit C
古典的情報 (2bit)
qubit C
EPR source
最大エンタングルメント でないとき忠実度の期待値
F (Fidelity)
• Alice , Bobが準備したエンタングルメント 状態の 2qubit
| Φ i ab = √
α| 0 i a | 1 i b + √
β| 1 i a | 0 i b ( α + β = 1 )
• Aliceが Bob
に送りたい状態| φ i c = c 0 | 0 i c + c 1 | 1 i c
• Alice
からの古典情報によってBob
が作り出した状態| φ
0i i
忠実度の期待値はF =
∑ 4 i = 1
P i |h φ
0i | φ i c | 2 ( P i
:状態| φ
0i i
を得る確率)= | c 0 | 4 + | c 1 | 4 + 4 | c 0 | 2 | c 1 | 2 √ αβ
= 1 − 2 | c 0 | 2 | c 1 | 2 (1 − C) (
C = 2 √ αβ )
α = β = 1 2
のとき(EPR pair) C = 1 = ⇒ F = 1
α , β
のときC < 1 = ⇒ F < 1
Concurrence
の導入Concurrence
とは・・・エンタングルメント の指標を表す。
Concurrence
の定義2qubit
の一般の状態|Φi
のConcurrence C( Φ )
は次の表し 方がある。•
任意の|Φi
に対して、Schmidt
の分解を行う。|Φi = √
α| 0
0i a | 0
0i b + √
β| 1
0i a | 1
0i b
( { | 0
0i a , | 1
0i a } , {| 0
0i b , | 1
0i b }
は正規直交基底)
→ C( Φ ) = 2 √ αβ
• θ = σ y ⊗ σ y K
とすると、( σ y
:y
方向のパウリ行列K
:複素共役をとる演算子)
→ C( Φ ) = |h Φ |θ| Φ i|
純粋状態と混合状態、密度演算子の導入
今までのように状態
| Φ i
にあるとき、量子系は純粋状態という。| Φ i = √
α| 00 i + √
β| 11 i
:純粋状態確率
q 1
で純粋状態| φ 1 i
に、確率q 2
で純粋状態| φ 2 i
にあるとき、量子系は混合状態と いう。ある
2qubit
で上記の混合状態での量子系における密度演算子ρ
は、ρ = q 1 | φ 1 ih φ 1 | + q 2 | φ 2 ih φ 2 |
一般に、確率
q i
で状態| φ i i
にあるような混合状態の密度演算子はρ = ∑
i
q i | φ i ih φ i |
純粋状態
( q i = 1
、状態は| φ i )
のときρ = | φ ih φ |
混合状態の
Concurrence
それぞれの状態の
Concurrence
の平均を混合状態のConcurrence
と仮定する。確率
q 1 = 1 2
で| Φ 1 i = √ 1
2 ( | 00 i + | 11 i )
確率q 2 = 1 2
で| Φ 2 i = √ 1
2 ( | 00 i − | 11 i ) ρ = 1
2 | Φ 1 ih Φ 1 | + 1
2 | Φ 2 ih Φ 2 | (1)
= 1
2 | 00 ih 00 | + 1
2 | 11 ih 11 | (2) (1)
の密度演算子のときC ( Φ 1 ) = 1
、C ( Φ 2 ) = 1
C ( ρ ) = q 1 C ( Φ 1 ) + q 2 C ( Φ 2 )
= 1
(2)
の密度演算子のときC (00) = 0
、C (11) = 0
C ( ρ ) = q 1 C (00) + q 2 C (11)
= 0
それぞれの結果からConcurrence
の値に矛盾が生じる。2
個のqubitがある混合状態にあるとき、エンタングルメント の指標である Concurrence
はどう定義するべきか?密度演算子
ρ
で記述される混合状態のConcurrence
を次のように定義する。C ( ρ ) = min ∑
i
q i C ( φ i )
ρ = ∑
i
q i | φ i ih φ i |
( min
はρ = ∑
q i | φ i ih φ i |
となるすべてのq i
と| φ i i
に関するものである)
この
min
を計算するのは容易ではないが、次の公式が成立することが分かっている。以降、次の公式を混合状態の
Concurrence
の定義として用いる。C ( ρ ) = max (
0 , √
λ 1 − √
λ 2 − √
λ 3 − √ λ 4
) ζ = ρθρθ = ρ (
σ y ⊗ σ y )
ρ ∗ (
σ y ⊗ σ y )
ただし 、
λ i
は演算子ζ
の固有値を大きい方から並べたものである。(λ 1 ≥ λ 2 ≥ λ 3 ≥ λ 4 )
2
準位原子の自然放射(温度 T = 0)
一つの2
準位原子を考える。時刻
t = 0
で励起状態にある原子は時刻t
にはある確率で自然放射し 、基底状態になる。他方、基底状態にある原子はそのまま基底状態に留まるとする。
密度演算子で書けば 、
ρ (0) = | 1 ih 1 |
↓
ρ ( t ) = (1 − e −Γ t ) | 0 ih 0 | + e −Γ t | 1 ih 1 | ρ (0) = | 0 ih 0 |
↓
ρ ( t ) = | 0 ih 0 |
ここで、
Γ 1
は励起状態の寿命である。一般に
2
準位原子の自然放射は次式で表される。ρ ( t ) = K 1 ( t ) ρ (0) K †
1 ( t ) + K 2 ( t ) ρ (0)K †
2 ( t ) = ∑
a = 1 , 2
K a ( t) ρ (0)K † a ( t )
K 1 ( t ) = | 0 ih 0 | + e
−Γ2 t| 1 ih 1 | K 2 ( t ) = √
1 − e −Γ t | 0 ih 1 |
K †
1 ( t ) = | 0 ih 0 | + e
−Γ2 t| 1 ih 1 | K †
2 ( t ) = √
1 − e −Γ t | 1 ih 0 |
K 1 ( t ) , K 2 ( t )
はクラウス演算子と呼ばれる。エンタングルメント の時間的変化
二つの
2
準位原子を考える。これらの間には相互作用はなく、独立に自然放射を行うと する。この時ρ ( t)
は次式で表される。ρ ( t) =
∑ 2 a , b = 1
K a ( t) ⊗ K b ( t ) ρ (0) K † a ( t ) ⊗ K †
b ( t )
相互作用なし 遠く離れている
⇐⇒
エンタングルメント の消滅および突然死
時刻
t = 0
で純粋状態、あるいは混合状態おける2
つのqubit
の時刻t
でのconcurrence
を調べる。純粋状態
ρ (0) = | 00 ih 00 | = ⇒ C( ρ ) = 0 ρ (0) = | 11 ih 11 | = ⇒ C( ρ ) = 0 ρ (0) = | ψ ih ψ | = ⇒ C( ρ ) = e −Γ t
(
| ψ i = √ 1
2
( | 01 i + | 10 i ) )
ρ (0) = | ψ 0 ih ψ 0 | = ⇒ C( ρ ) = e − 2 Γ t
(
| ψ 0 i = √ 1
2
( | 00 i + | 11 i )
)
混合状態
ρ (0) = (1 − a) | 00 ih 00 | + a | ψ ih ψ |
(
| ψ i = √ 1
2
( | 01 i + | 10 i ) )
= ⇒ C( ρ ) = max {
ae −Γ t , 0 }
(0 ≤ a ≤ 1)
0 0.2 0.4 0.6 0.8 1
0 0.5 1 1.5 2 2.5 3
Concurrence
time
a=1 a=0
ρ (0) = (1 − a) | 11 ih 11 | + a | ψ ih ψ |
(
| ψ i = √ 1
2
( | 01 i + | 10 i ) )
= ⇒ C( ρ ) = max {
ae −Γ t − 2e −Γ t √
(1 − a)(1 − e −Γ t )(ae −Γ t − e −Γ t + 1) , 0 }
0 0.2 0.4 0.6 0.8 1
0 0.5 1 1.5 2 2.5 3
Concurrence
time
a=1 a=0
0 0.2 0.4 0.6 0.8 1
0.5 0 1.5 1
2.5 2 3
0 0.2 0.4 0.6 0.8 1
a time
C
上記の混合状態の結果から、
2
原子が励起した状態| 11 i
と、エンタングルした状態| ψ i
との混合状態で はConcurrence
の値は有限のt
で0
になることがある。= ⇒
突然死!2
原子が基底状態にある状態| 00 i
と、エンタングルした状態| ψ i
との混合 状態ではConcurrence
の値は指数関数的に0
に収束する。(
| ψ i = √ 1
2
( | 01 i + | 10 i ) )
「エンタングルメント の突然死」は
J.H.Eberly
とTing Yu
によって発見された。(2006
年)(次ページ参照)
•
上記の混合状態の組み合わせは私達が考えたものだが、同様に突然死が起きる場合 も確認できた。J.H.Eberly
とTing Yu
の計算ρ (0) = 1 3 (a | 00 ih 00 |
+ | 01 ih 01 | + | 01 ih 10 | + | 10 ih 01 | + | 10 ih 10 | + (1 − a) | 11 ih 11 | )
C( ρ ) = max { 2
3 e −Γ t { 1 − √
(1 − a) { (1 − a)(1 − e −Γ t ) 2 + 2(1 − e −Γ t ) + a }} , 0
}
a =1
C( ρ ) = max { 2
3 e −Γ t , 0 }
a =0
C( ρ ) = max { 2
3 e −Γ t { 1 − √
(1 − e −Γ t ) 2 + 2(1 − e −Γ t ) } , 0 }
Γ t = log (
1 + √ 1
2
)
で突然死する。
0 0.2 0.4 0.6 0.8 1
0.5 0 1.5 1
2.5 2 3
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
a
timeC
a=1 a=0
有限温度でのエンタングルメント
有限温度では一般に
2
準位原子の光の放射と吸収は次式で表される。ρ ( t ) = ∑
a = 1 , 2 , 3 , 4
K a ( t ) ρ (0)K † a ( t )
K 1 ( t ) = √
1 − P( | 0 ih 0 | + e
−Γ2 t| 1 ih 1 | ) K 3 ( t ) = √
P(e
−Γ2t| 0 ih 0 | + | 1 ih 1 | ) K 2 ( t ) = √
1 − P( √
1 − e −Γ t | 0 ih 1 | ) K 4 ( t ) = √
P( √
1 − e −Γ t | 1 ih 0 | )
上式を用いて
ρ (0) = | 0 ih 0 |, ρ (0) = | 1 ih 1 |
の時刻t → ∞
での密度演算子を求めると、ρ ( ∞ ) = (1 − P) | 0 ih 0 | + P | 1 ih 1 |
となる。 ここで
P
は、P = P 1 = e
−1 kT
e
−0 kT
+ e
−1 kT
= e
−(1−0) kT
1 + e
−(1−0) kT
( P 1 : t → ∞
で励起状態に留まる確率。)
T = 0 → ρ ( ∞ ) = | 0 ih 0 |
T = ∞ → ρ ( ∞ ) = 1 2 ( | 0 ih 0 | + | 1 ih 1 | )
エンタングルメント の
Concurrence (有限温度)
有限温度では一般に二つの
2
準位原子の光の放射と吸収は次式で表される。ρ ( t ) =
∑ 4 a , b = 1
K a ( t ) ⊗ K b ( t ) ρ (0)K † a ( t ) ⊗ K †
b ( t)
K 1 ( t ) = √
1 − P( | 0 ih 0 | + e
−Γ2 t| 1 ih 1 | ) K 3 ( t ) = √
P(e
−Γ2t| 0 ih 0 | + | 1 ih 1 | ) K 2 ( t ) = √
1 − P( √
1 − e −Γ t | 0 ih 1 | ) K 4 ( t ) = √
P( √
1 − e −Γ t | 1 ih 0 | )
時刻
t = 0
で純粋状態ρ (0) = | ψ ih ψ |
(
| ψ i = √ 1
2
( | 01 i + | 10 i ) )
= ⇒ C( ρ ) = max {
e −Γ t − 2 P(1 − P)(1 − e −Γ t )
√
1 + e − 2 Γ t + e −Γ t ( 1
P(1 − P) − 2) , 0 }
T = 0
C( ρ ) = e −Γ t T = ∞
C( ρ ) = 1 2 (e −Γ t + 1) 2 − 1 Γ t = log
(
√ 1 2 − 1
)
で突然死する。 0
0.2 0.4 0.6 0.8 1
0 0.5 1 1.5 2 2.5 3
Concurrence
time
P=1/2 P=1/3 P=1/6 P=0
結論
•
量子エンタングルメント を利用することで、様々な量子情報的な作業を行うことが 出来る。(量子テレポーテーションや Super dense coding
など)
•
量子テレポーテーションの忠実度F
はエンタングルメント の指標であるConcur- rence
と同じになる。•
温度T = 0
で独立に自然放射を行う2
準位原子を考えたとき、励起した状態
| 11 i
とエンタングルした状態| ψ i
との混合状態ではConcurrence
の 値は有限のt
で0
になることがあり、突然死する。(
| ψ i = √ 1
2
( | 01 i + | 10 i ) )
•
有限温度では時刻t = 0
で純粋状態でもConcurrence
の値は有限のt
で0
になり、突然死する。
補足1
(量子テレポーテーション )
•
ベル基底
| Φ 1 i = √ 1 2
( | 00 i + | 11 i )
| Φ 2 i = √ 1
2 ( | 00 i − | 11 i )
| Φ 3 i = √ 1
2 ( | 01 i + | 10 i )
| Φ 4 i = √ 1
2 ( | 10 i − | 10 i )
•
ユニタリー変換
U 1 = 1
U 2 = | 0 ih 0 | − | 1 ih 1 | U 3 = | 1 ih 0 | + | 0 ih 1 | U 4 = U 2 U 3
•
量子テレポーテーション| Φ i abc = | Φ 1 i ac
( C
0√
√ α
2 | 0 i b + C
1√ β
√ 2 | 1 i b
)
+ | Φ 2 i ac
( C
0√
√ α
2 | 0 i b − C
1√ β
√ 2 | 1 i b
) + | Φ 3 i ac
( C
0√
√ β
2 | 1 i b + C
1√ √ α
2 | 0 i b
)
+ | Φ 4 i ac
( C
0√
√ β
2 | 1 i b − C
1√ √ α
2 | 0 i b
)
Aliceが Bell
測定を行う2bit(4
通り)のいづれかの測定結果を得る→その結果に応じて
Bob
はユニタリー変換を行う。補足2
(Concurrence)
• Schmidt
の標準形を用いる場合| Φ i = α 00 | 00 i + α 01 | 01 i + α 10 | 10 i + α 11 | 11 i
= ( α 00 | 0 i a + α 10 | 1 i a ) | 0 i b + ( α 01 | 0 i a + α 11 | 1 i a ) | 1 i b
= √
|α 00 | 2 + |α 10 | 2
α
00√ | |α 0 i
a+α
10| 1 i
a00
|
2+|α
10|
2
| 0 i b + √
|α 01 | 2 + |α 11 | 2
α
01√ | |α 0 i
a+α
11| 1 i
a01
|
2+|α
11|
2
| 1 i b
= √
|α 00 | 2 + |α 10 | 2 | 0
0i a | 0
0i b + √
|α 01 | 2 + |α 11 | 2 | 1
0i a | 1
0i b
= √
α| 0
0i a | 0
0i b + √
β| 1
0i a | 1
0i b
C( Φ ) = 2 √ αβ
• y
方向のパウリ行列を用いる場合θ = σ y ⊗ σ y K
= ( −| 00 ih 11 | + | 01 ih 10 | + | 10 ih 01 | − | 11 ih 00 | ) K
| Φ i = α 00 | 00 i + α 01 | 01 i + α 10 | 10 i + α 11 | 11 i
= √
α| 0
0i a | 0
0i b + √
β| 1
0i a | 1
0i b
C( Φ ) = |h Φ |θ| Φ i| = − 2 √
αβ = 2 √
αβ
補足
3(密度演算子)
ある状態
| Φ i
が次のように表されている.| Φ i = c 1 | 1 i + c 2 | 2 i + · · · + c a | a i + · · · =
∑ n a = 1
c a | a i
正規直交基底で測定をすると、結果a
を得る確率はP a = | c a | 2 = |h a | Φ i| 2
確率
q 1
で純粋状態| φ i 1
、確率q 2
で純粋状態| φ i 2
のような混合状態にある量子系を正規 直交基底で測定を行い、結果a
を得る確率P a = q 1 |h a | φ 1 i| 2 + q 2 |h a | φ 2 i| 2
= q 1 h a | φ 1 ih φ 1 | a i + q 2 h a | φ 2 ih φ 2 | a i
= h a | (q 1 | φ 1 ih φ 1 | + q 2 | φ 2 ih φ 2 | ) | a i
= h a |ρ| a i
ρ ≡ q 1 | φ 1 ih φ 1 | + q 2 | φ 2 ih φ 2 |
確率